The Influence of Previous Mountain Pine Beetle (Dendroctonus ponderosae) Activity on the 1988 Yellowstone Fires

We examined the historical record of mountain pine beetle (Dendroctonus ponderosae Hopkins) activity within Yellowstone National Park, Wyoming, for the 25-years period leading up to the 1988 Yellowstone fires (1963–86) to determine how prior beetle activity and the resulting tree mortality affected the spatial pattern of the 1988 Yellowstone fires. To obtain accurate estimates of our model parameters, we used a Markov chain Monte Carlo method to account for the high degree of spatial autocorrelation inherent to forest fires. Our final model included three statistically significant variables: drought, aspect, and sustained mountain pine beetle activity in the period 1972–75. Of the two major mountain pine beetle outbreaks that preceded the 1988 fires, the earlier outbreak (1972–75) was significantly correlated with the burn pattern, whereas the more recent one (1980–83) was not. Although regional drought and high winds were responsible for the large scale of this event, the analysis indicates that mountain pine beetle activity in the mid-1970s increased the odds of burning in 1988 by 11% over unaffected areas. Although relatively small in magnitude, this effect, combined with the effects of aspect and spatial variation in drought, had a dramatic impact on the spatial pattern of burned and unburned areas in 1988.     

Physiological model using diffuse reflectance spectroscopy for nonmelanoma skin cancer diagnosis

Diffuse reflectance spectroscopy (DRS) is a noninvasive, fast, and low‐cost technology with potential to assist cancer diagnosis. The goal of this study was to test the capability of our physiological model, a computational Monte Carlo lookup table inverse model, for nonmelanoma skin cancer diagnosis. We applied this model on a clinical DRS dataset to extract scattering parameters, blood volume fraction, oxygen saturation and vessel radius. We found that the model was able to capture physiological information relevant to skin cancer. We used the extracted parameters to classify (basal cell carcinoma [BCC], squamous cell carcinoma [SCC]) vs actinic keratosis (AK) and (BCC, SCC, AK) vs normal. The area under the receiver operating characteristic curve achieved by the classifiers trained on the parameters extracted using the physiological model is comparable to that of classifiers trained on features extracted via Principal Component Analysis. Our findings suggest that DRS can reveal physiologic characteristics of skin and this physiologic model offers greater flexibility for diagnosing skin cancer than a pure statistical analysis. Physiological parameters extracted from diffuse reflectance spectra data for nonmelanoma skin cancer diagnosis.     

Effect of anode angle on photon beam spectra and depth dose characteristics for X-RAD320 orthovoltage unit

Background

In radiation therapy with orthovoltage units, the tube design has a crucial effect on its dosimetric features.

Aim

In this study, the effect of anode angle on photon beam spectra, depth dose and photon fluence per initial electron was studied for a commercial orthovoltage unit of X-RAD320 biological irradiator.

Materials and methods

The MCNPX MC code was used for modeling in the current study. We used the Monte Carlo method to model the X-RAD320 X-ray unit based on the manufacturer provided information. The MC model was validated by comparing the MC calculated photon beam spectra with the results of SpekCalc software. The photon beam spectra were calculated for anode angles from 15 to 35 degrees. We also calculated the percentage depth doses for some angles to verify the impact of anode angle on depth dose. Additionally, the heel effect and its relation with anode angle were studied for X-RAD320 irradiator.

Results

Our results showed that the photon beam spectra and their mean energy are changed significantly with anode angle and the optimum anode angle of 30 degrees was selected based on less heel effect and appropriate depth dose and photon fluence per initial electron.

Conclusion

It can be concluded that the anode angle of 30 degrees for X-RAD320 unit used by manufacturer has been selected properly considering the heel effect and dosimetric properties.     

A theoretical model for the association of amphiphilic transmembrane peptides in lipid bilayers

A theoretical model is proposed for the association of trans-bilayer peptides in lipid bilayers. The model is based on a lattice model for the pure lipid bilayer, which accounts accurately for the most important conformational states of the lipids and their mutual interactions and statistics. Within the lattice formulation the bilayer is formed by two independent monolayers, each represented by a triangular lattice, on which sites the lipid chains are arrayed. The peptides are represented by regular objects, with no internal flexibility, and with a projected area on the bilayer plane corresponding to a hexagon with seven lattice sites. In addition, it is assumed that each peptide surface at the interface with the lipid chains is partially hydrophilic, and therefore interacts with the surrounding lipid matrix via selective anisotropic forces. The peptides would therefore assemble in order to shield their hydrophilic residues from the hydrophobic surroundings. The model describes the self-association of peptides in lipid bilayers via lateral and rotational diffusion, anisotropic lipid-peptide interactions, and peptide-peptide interactions involving the peptide hydrophilic regions. The intent of this model study is to analyse the conditions under which the association of trans-bilayer and partially hydrophilic peptides (or their dispersion in the lipid matrix) is lipid-mediated, and to what extent it is induced by direct interactions between the hydrophilic regions of the peptides. The model properties are calculated by a Monte Carlo computer simulation technique within the canonical ensemble. The results from the model study indicate that direct interactions between the hydrophilic regions of the peptides are necessary to induce peptide association in the lipid bilayer in the fluid phase. Furthermore, peptides within each aggregate are oriented in such a way as to shield their hydrophilic regions from the hydrophobic environment. The average number of peptides present in the aggregates formed depends on the degree of mismatch between the peptide hydrophobic length and the lipid bilayer hydrophobic thickness: The lower the degree of mismatch is the higher this number is. Received: 30 December 1996 / Accepted: 9 May 1997     

Economic analysis of royalactin production under uncertainty: Evaluating the effect of parameter optimization

Royalactin is a protein with several different potential uses in humans. Research, in insects and in mammalian cells, has shown that it can accelerate cell division and prevent apoptosis. The method of action is through the use of the epidermal growth factor receptor, which is present in humans. Potential use in humans could be to lower cholesterolemic levels in blood, and to elicit similar effects to those seen in bees, e.g., increased lifespan. Mass production of Royalactin has not been accomplished, though a recent article presented a Pichia pastoris fermentation and recovery by aqueous two‐phase systems at laboratory scale as a possible basis for production. Economic modelling is a useful tool with which compare possible outcomes for the production of such a molecule and in particular, to locate areas where additional research is needed and optimization may be required. This study uses the BioSolve software to perform an economic analysis on the scale‐up of the putative process for Royalactin. The key parameters affecting the cost of production were located via a sensitivity analysis and then evaluated by Monte Carlo analysis. Results show that if titer is not optimized the strategy to maintain a low cost of goods is process oriented. After optimization of this parameter the strategy changes to a product‐oriented and the target output becomes the critical parameter determining the cost of goods. This study serves to provide a framework for the evaluation of strategies for future production of Royalactin, by analyzing the factors that influence its cost of manufacture. © 2015 American Institute of Chemical Engineers Biotechnol. Prog., 31:744–749, 2015     

Health Risk Assessments Based on Existing Data of Arsenic,Chromium, Lead,and Zinc in China's Air

Health risk assessments of As, Cr, Pb, and Zn in air in China were carried out from a national perspective. Concentrations of the four metals in air were gathered from published literature. Health risk assessment models recommended by the U.S. Environmental Protection Agency (USEPA) were applied to quantitatively characterize risk values of the metals by considering different age groups. Results showed that the atmospheric metal concentrations in most regions of China were lower than their inhalation reference concentrations provided by the USEPA, or limited values provided by China. Non-carcinogenic hazard values of the four metals in different regions of China were all lower than 1, indicating no or low non-carcinogenic hazard. However, the probabilities of carcinogenic risk values for As and Cr exceeding the acceptable value (1.0E-05) were 9% and 10% on a national scale, respectively. Some regions, such as Shanxi, Chongqing, Liaoning, and Shanghai, should be studied in more detail. Despite uncertainties, the results of this study provide information on the health risk of As, Cr, Pb, and Zn in air in China.     

Monte Carlo simulation of the dose distribution of ICRP adult reference computational phantoms for acquisitions with a 320 detector-row cone-beam CT scanner

PurposeThe purpose of this study was to develop and validate a Monte Carlo (MC) simulation tool for patient dose assessment for a 320 detector-row CT scanner, based on the recommendations of International Commission on Radiological Protection (ICRP). Additionally, the simulation was applied on four clinical acquisition protocols, with and without automatic tube current modulation (TCM).MethodsThe MC simulation was based on EGS4 code and was developed specifically for a 320 detector-row cone-beam CT scanner. The ICRP adult reference phantoms were used as patient models. Dose measurements were performed free-in-air and also in four CTDI phantoms: 150 mm and 350 mm long CT head and CT body phantoms. The MC program was validated by comparing simulations results with these actual measurements acquired under the same conditions. The measurements agreed with the simulations across all conditions within 5%. Patient dose assessment was performed for four clinical axial acquisitions using the ICRP adult reference phantoms, one of them using TCM.ResultsThe results were nearly always lower than those obtained from other dose calculator tools or published in other studies, which were obtained using mathematical phantoms in different CT systems. For the protocol with TCM organ doses were reduced by between 28 and 36%, compared to the results obtained using a fixed mA value.ConclusionsThe developed simulation program provides a useful tool for assessing doses in a 320 detector-row cone-beam CT scanner using ICRP adult reference computational phantoms and is ready to be applied to more complex protocols.     

Oxygenation measurement by multi‐wavelength oxygen‐dependent phosphorescence and delayed fluorescence: catchment depth and application in intact heart

Oxygen delivery and metabolism represent key factors for organ function in health and disease. We describe the optical key characteristics of a technique to comprehensively measure oxygen tension (PO₂) in myocardium, using oxygen‐dependent quenching of phosphorescence and delayed fluorescence of porphyrins, by means of Monte Carlo simulations and ex vivo experiments. Oxyphor G2 (microvascular PO₂) was excited at 442 nm and 632 nm and protoporphyrin IX (mitochondrial PO₂) at 510 nm. This resulted in catchment depths of 161 (86) µm, 350 (307) µm and 262 (255) µm respectively, as estimated by Monte Carlo simulations and ex vivo experiments (brackets). The feasibility to detect changes in oxygenation within separate anatomical compartments is demonstrated in rat heart in vivo.

Schematic of ex vivo measurements.     

Replica-permutation method to enhance sampling efficiency

In order to predict the native structures of proteins and peptides and to investigate the functions of these biomolecules, it is essential to realise efficient sampling in the conformational space. We had recently proposed a new generalised-ensemble algorithm, which is referred to as the replica-permutation method (RPM), to sample the conformational space efficiently. We introduce this RPM and demonstrate its usefulness by applying to three systems: particles in a double-well potential energy, Met-enkephalin in a vacuum, and a C-peptide analogue of ribonuclease A in explicit water. Replica-exchange simulations were performed to compare their results with the results of the replica-permutation simulations. It is shown that the RPM sampled not only the temperature space but also the conformational space more efficiently than the replica-exchange method. The folding pathway of C-peptide is also presented.     

Cation redistribution upon dehydration of Na58Y faujasite zeolite: a joint neutron diffraction and molecular simulation study

Using neutron scattering and Monte Carlo simulation, we investigate the distribution of cations in Na₅₈Y faujasite upon (de)hydration. We introduce a new method for the assignment of cations to specific sites in molecular simulations from their local environment. This allows us to bypass the need of the coordinates of crystallographic sites, which vary as water adsorption induces changes in the zeolite framework structure. Although the agreement between experiments and simulation is excellent at high temperature, some differences are observed below 150°C. We show that these differences are due to the presence of water and that temperature itself as well as adsorption-induced deformation of the framework play a less important role. We demonstrate the migration of sodium to sites III upon water adsorption, not observed for other Si:Al ratios.